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Numerous proteins require cofactors to be active. Computer simulations suggest that 
cooperative interaction networks achieve optimal cofactor binding. There is a need for the 
experimental identification of the residues crucial for stabilizing these networks and thus for 
cofactor binding. Here we investigate the electron transporter flavodoxin, which contains 
flavin mononucleotide as non-covalently bound cofactor. We show that after binding flavin 
mononucleotide with nanomolar affinity, the protein relaxes extremely slowly (time constant 
~5 days) to an energetically more favourable state with picomolar-binding affinity. Rare 
small-scale openings of this state are revealed through H/D exchange of N(3)H of flavin. 
We find that H/D exchange can pinpoint amino acids that cause tight cofactor binding. These 
hitherto unknown residues are dispersed throughout the structure, and many are located 
distantly from the flavin and seem irrelevant to flavodoxin's function. Quantification of the 
thermodynamics of ligand binding is important for understanding, engineering, designing and 
evolving ligand-binding proteins. 
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Many proteins bind ligands, which usually contact a rather 
limited number of amino -acid residues. However, their 
binding effects can propagate to residues that have no 
direct interaction with the ligand and as a result enable biologi- 
cal phenomena such as allosteric regulation and signal transduc- 
tion^'^. During binding, a ligand selects the protein conformers that 
best match its structural and dynamic properties among the pool 
of conformers intrinsically accessible to the protein in unliganded 
form^. Cofactors form a special category of ligands. A cofactor is a 
non-protein chemical compound that is bound to a protein and is 
required for the proteins biological activity. On binding, the prop- 
erties of both cofactor and protein can mutually alter, as happens, 
for instance, in case of hemo- or flavoproteins, thereby facilitating 
important life processes like electron transfer. 

Dynamic fluctuations are fundamentally linked to protein func- 
tion. Fluorescence and nuclear magnetic resonance (NMR) spec- 
troscopy have illuminated fast internal protein dynamics and their 
biological purpose. For example, NMR relaxation methodologies, 
which reveal protein dynamics on the millisecond timescale and 
below, enable characterization of dynamic features of enzymes rele- 
vant for catalysis^'^. On ligand binding, significant changes in fast 
internal dynamics can occur. This perturbation might be local, or 
transmitted to sites remote from the binding site^. Usually, binding 
of a ligand reduces protein flexibility. However, nature sometimes 
employs increases in the residual entropy of proteins, manifested 
in their conformational dynamics, as a thermodynamic mechanism 
for allosteric responses to ligand binding. Hence, enhanced internal 
protein dynamics on ligand binding have been observed as welF'^. 

Computer simulations suggest that propagation of ligand-bind- 
ing interactions and protein stability are linked. This link is medi- 
ated through distribution of stability within different structural 
elements of a protein, and interactions existing between them^. By 
means of these interaction networks, ligand-induced stabilization 
can transmit throughout a protein. Although considerable progress 
has been made in computational studies of ligand binding^^"^^, 
a strong need exists for experimental data that probe and verify, 
at the residue level, the effects ligand binding have on protein 
thermodynamic stability. 

Flavoproteins are cofactor- containing proteins that participate 
in various biological functions, such as photosynthesis, cellular res- 
piration, DNA repair, apoptosis, and nerve signal transmission^^. 
Flavodoxins are monomeric electron transfer proteins and contain 
flavin mononucleotide (FMN) that can exist as oxidized quinone, 
one electron-reduced semiquinone and two-electron reduced 
hydroquinone^^'^^. 

The redox potentials of flavodoxins are affected by several 
factors, including a conserved tyrosine that tightly packs onto the 
isoalloxazine ring of FMN and differentially interacts with it on 
changing redox states, stabilization of the semiquinone relative to 
the oxidized form (because it forms a new hydrogen bond with a 
main chain group), repelling interactions of the negatively charged 
isoalloxazine ring of two-electron-reduced FMN with both acidic 
residues on the periphery of the protein and the negatively charged 
phosphate of FMN, and unfavourable burial inside the hydrophobic 
interior of flavodoxin of this uncompensated negative charge^^"^^. 
As a result, flavodoxin can transfer electrons between various physi- 
ological redox partners. Strong binding of FMN is essential for fla- 
vodoxins to be functionally active, and these proteins are remarkably 
rigid on the ns-timescale and below^^"^^. The isoalloxazine moiety 
and 5'-phosphate group of FMN contribute almost equally to the 
free energy associated with cofactor binding, whereas the ribityl 
has barely any impact^^. Currently, amino-acid residues crucial for 
stabilizing FMN-protein interactions are unknown. Thus, whereas 
the redox potential-affecting contributions of flavodoxin have been 
revealed, the structural determinants that govern strong cofactor 
binding are still obscure. 



Use of H/D exchange detected by NMR spectroscopy can iden- 
tify amino-acid residues that are required for tight cofactor binding. 
Because of an unfolding process, a residue can arrive in an open 
or exchange -competent form from which amide proton exchange 
with deuterons takes place when a protein is dissolved in deuterium 
oxide. For some residues, the dominant process that causes amide 
exchange to occur is a highly local opening of the native protein 
structure. For other residues, exchange only happens if rare open- 
ing events like subglobal or global protein unfolding occur. These 
unfolding processes typically occur on timescales ranging from 
seconds to many days^^'^^. Using H/D exchange data, thermody- 
namic stabilities against unfolding can be quantified, leading to first 
efforts to address the influences ligand binding has on local protein 
stabilities (see example in refs 26-29). 

To correctly interpret H/D exchange data of a cofactor- containing 
protein, knowledge about the (un) folding and local thermodynamic 
stabilities of the corresponding apoprotein is required, as well as a 
thorough understanding of the role cofactor binding has in pro- 
tein folding. We demonstrated that the in vitro folding of apoflavo- 
doxin (that is, flavodoxin without FMN) occurs autonomously and 
spontaneously and involves two intermediates^^"^^. Non-covalent 
binding of FMN to native apoflavodoxin is the last step in flavo- 
doxin folding, and during global unfolding of flavodoxin, release of 
FMN happens first^^ NMR data showed that native apoflavodoxin 
strongly resembles flavodoxin, except for dynamic disorder in the 
flavin-binding region^^. 

Here we report the effects cofactor binding has on local ther- 
modynamic stabilities of a flavodoxin from Azotobacter vinelandii. 
We demonstrate that after reconstitution from apoprotein and 
FMN, flavodoxin relaxes extremely slowly to an energetically more 
favourable state with two orders of magnitude better FMN binding. 
We apply H/D exchange to investigate local thermodynamic sta- 
bilities of relaxed flavodoxin, which has picomolar binding affinity 
for FMN. We show that N(3)H of flavin is an excellent probe to fol- 
low opening events in flavin-binding sites and identify amino-acid 
residues that cause tight binding of the cofactor. 

Results 

Determination of the strength of FMN binding. We determined 
the equilibrium constant for dissociation of FMN from flavodoxin 
(that is, Kiy) according to a well-established procedure in which 
aliquots of apoprotein are titrated to an FMN solution (Supple- 
mentary Methods). Use is made of the severe quenching of the 
fluorescence of FMN on its binding to apoprotein^ Quenching 
happens within seconds^ ^ After each addition of apoprotein, the 
system was allowed to eqmlibrate for 5min in the dark before 
recording of FMN fluorescence (Fig. la). We derive that is 
(3.51±0.20)-10"i^M, which is fully consistent with Kp-values 
reported for other flavodoxins^^. 

To quantify the rate constant for FMN release, /Coff, we followed 
the increase in flavin fluorescence due to cofactor dissociation after 
dilution of flavodoxin in buffer. Reestablishment of equilibrium 
between apo- and holo-protein proceeds mono -exponentially as 
function of time, allowing determination of /Coff, as well as (Sup- 
plementary Methods) Figure lb reports four relaxation traces as 
well as the global fit of equation (SI 1) of Supplementary Methods to 
the data. We thus derive that k^^^is (1.05±0.04)-10"6s" ^ Remark- 
ably now Ky) is established to be (3.82±0.12)-10" ^^M; that is, 
FMN binds two orders of magnitude better to apoprotein than the 
titration procedure reports. 

Fluorescence spectroscopy is used to elucidate the source for the 
difference in observed. Fluorescence of 5 \\M flavodoxin shows 
a peak at 525 nm characteristic for free FMN (Fig. 2a), which is 
dissociated from holoprotein as dictates. On adding excess 
apoprotein, free FMN is no longer detected. In case of 5 |lM recon- 
stituted flavodoxin, obtained by mixing apoflavodoxin and FMN 
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Figure 1 1 Equilibrium constant for dissociation of FMN from flavodoxin 
and corresponding rate constant, (a) Determination of /Cp using the 
quenching of FMN fluorescence on binding of this cofactor to apoprotein 
at 25 °C. A 188-nM FMN solution is titrated with aliquots of a 4.4-|iM 
apoflavodoxin solution. One unit of fluorescence corresponds to 
0.22 nM FMN. Equation (S8) of Supplementary Methods is fitted to 
the resulting fluorescence intensity data and /Cq is determined to be 
(3.51 ±0.20)-10~^^ M. The two linear components of this equation are 
shown as dashed lines. The inset highlights the curved part of the titration 
data, (b) On dilution of flavodoxin in buffer at 25 °C, FMN is released 
and fluorescence increases accordingly (fluorescence emission of FMN 
is followed at 525 nm, while excitation occurred at 450 nm). One unit 
of fluorescence corresponds to 5.0 nM FMN. After dilution of 5 (open 
circles), 10 (closed circles), 20 (open squares) and 50 (closed squares) 
jil of 362jiM flavodoxin in 2,000 |il buffer, respectively, reestablishment 
of the equilibrium between apo- and holoprotein proceeds exponentially 
as function of time. Equation (S11) of Supplementary Methods is 
globally fitted to the four relaxation traces and shows that /Cgff is 
(1.05±0.04)-10-6s-^ and is (3.82±0.12)-10-^2 



at equimolar ratio higher FMN fluorescence is observed. On 
adding excess apoprotein a peak is detected at 512nm (Fig. 2a). 
This 13nm blue-shifted fluorescence cannot be due to free flavin, 
but must arise from FMN bound to reconstituted flavodoxin. 
Clearly, reconstituted protein differs from flavodoxin, explaining 
the difference in Kj^-values observed. Inspection of FMN fluores- 
cence of reconstituted flavodoxin reveals that disappearance of 
fluorescence proceeds extremely slowly with time on relaxation of 
holoprotein (Fig. 2b). 

When freshly reconstituted flavodoxin is subjected to a quick 
cycle of two-electron reduction and immediate re-oxidation, we 
detect full quenching of FMN fluorescence, which is the hallmark 
of relaxed flavodoxin. This observation reveals that tightening of 
cofactor binding can happen rapidly. On reduction of flavodoxin, 
N(5) of the flavin becomes protonated. This N(5)H proton forms 
a hydrogen bond with a specific carbonyl oxygen of the protein 
backbone, causing the corresponding peptide bond to flip^^. This 
redox-modulated conformational change is the most likely source 
that causes rapid tightening of flavin binding. This rapid tightening 
is of functional relevance, because flavodoxins act in the reducing 
environment of the cytoplasm as one -electron carriers by cycling 
between the semiquinone and hydroquinone forms of FMN. 

To identify flavodoxin residues involved in slow structural 
relaxation after cofactor binding, ^^N-labelled apoflavodoxin and 
^^N-labelled FMN were mixed. Subsequently, a series of ^H-^^N 
heteronuclear single -quantum coherence (HSQC) spectra was 
recorded, which reveals that 54 backbone amide cross-peaks of this 
179-residue protein shift during a period of 58 days (Fig. 3). 

Quantification of thermodynamic stabilities. To reveal how FMN 
locally affects the stabilities of flavodoxins residues against exposure 
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Figure 2 | Fluorescence emission of flavodoxin and of reconstituted 
holoprotein. (a) Spectrum-2 is of 5jiM flavodoxin in buffer, whereas 
spectrum-1 is of 5jil\/l reconstituted flavodoxin, which is obtained by 
mixing apoflavodoxin and FMN at equimolar ratio. On subsequent addition 
of excess apoflavodoxin to a final concentration of 375 nM, fluorescence 
emission of both samples changes differently (spectrum-4 and spectrum-3, 
respectively). This addition reveals that fluorescence of FMN in flavodoxin 
is fully quenched, whereas in freshly reconstituted holoprotein, it is 
not. (b) Time-dependent change in fluorescence emission of 5jiM 
reconstituted flavodoxin (excitation 450 nm, emission 514 nm). One 
unit of fluorescence corresponds to 0.88 nM FMN. Single-exponential 
fitting of the data shows that the relaxation time is 4.9 + 0.2 days. 



to solvent, we measured H/D exchange rates (that is, k^^) of back- 
bone amides of apo- as well as of non- reconstituted holo-protein 
(Supplementary Table SI). Lyophilized ^^N-labelled protein was 
dissolved into D2O containing 100 mM potassium pyrophosphate at 
25 °C, pD 6.3. Each sample was transferred immediately into a NMR 
instrument and subsequently a series of ^H-^^N HSQC spectra was 
recorded. To obtain k^^, we analysed time -dependent maximal cross- 
peak intensities of backbone amides (Supplementary Methods). 
Figure 4a- e show typical exchange curves for flavodoxin. 

Quantitative interpretation of H/D exchange is possible using a 
simple model in which an open or exchange-competent form and a 
closed or exchange-incompetent form of a protein at the site of a par- 
ticular amide proton interconvert (Supplementary Methods)^^. Two 
limiting situations exist in H/D exchange. In case of EXl exchange 
behaviour, the exchange rate k^^ does not depend on pD and equals 
the rate constant for local opening of the protein structure, k^^. In 
case of EX2 exchange behaviour, log(A:ex) depends linearly on pD 
with a slope of 1 and /c^x-values can be directly converted into free 
energy differences between open and closed states, that is, AGop-values 
(Fig. 5a,b; Supplementary Methods). EX2 exchange behaviour has 
been shown to apply for apoflavodoxin at pD 6.3 (ref. 33), and thus 
we can extract AGop-values from the H/D exchange rates of the 
backbone amides of apoflavodoxin (Fig. 6a; Supplementary Table 1). 
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Figure 3 | Slow structural relaxation within reconstituted flavodoxin. (a) Overlay of ^H-^^N HSQC spectra recorded 33 min (blue) and 58 days (red) 
after addition of ^^N-labelled FMN to ^^N-labelled apoflavodoxin. Flavodoxin concentration is 2.2 mM in 100 mM potassium pyrophosphate, 90% 
H20/10% D2O, pH 6.0, at 25 °C. Cross-peaks that shift prominently are highlighted by black boxes. The inset zooms in on the area that is outlined by 
a grey box. (b) Amino-acid residues of which the corresponding cross-peaks shift by 0.003 p. p.m. or more are coloured red in the cartoon model of 
flavodoxin. 



To determine the mechanism by which amides of flavodoxin 
exchange, we measured H/D exchange at pD 6.2, 6.3, 7.7 and 8.8 
(Supplementary Table 2, 25 °C). The strong pD-dependence of 
the exchange data collected shows that at pD 6.3, EX2 conditions 
apply (Supplementary Methods, and Fig. 5b as example). This 
observation implies that at this pD, one can directly convert amide 
exchange rates of flavodoxin into corresponding AGop -values 
(Fig. 6b; Supplementary Table SI). 

FMN-induced stabilization varies widely across holoprotein. The 

H/D exchange data set obtained of apo- and holo-protein (Fig. 4; 
Supplementary Tables SI and S2) allows assessment of the effects 
FMN binding has on protecting residues against exchange. The term 
AAGop = (AGop(flavodoxin)-AGop(apoflavodoxin)) quantifies the 
change in stability against exposure to solvent of the residues of apo- 
flavodoxin as a result of FMN binding. For 87 residues of flavodoxin, 
we obtained AAGop-values, because the corresponding AG^p-values 
of both apo- and holo-protein are available (Supplementary 
Table SI; the standard deviation of AGop is 0.35kcalmol~^ (Sup- 
plementary Methods), and consequently, the standard deviation of 
AAGop is 0.50kcalmori). 

Figure 7 shows that the magnitudes of AAGop-values vary 
widely across flavodoxin. These magnitudes do not correlate with 
AGop(apoflavodoxin) (Supplementary Fig. SI). For 30 residues (col- 
oured pale blue in Fig. 7), the corresponding AAGop-values equal 
within error 0 kcalmol"^ (that is, -0.5< AAGop < 0.5 kcalmol"^). 
These residues reside in parts of the protein where binding of 
FMN apparently does not affect protein dynamics relevant for H/D 
exchange. Many of these 30 residues are located within or close to 
the part of the parallel p-sheet that is opposite to the site where 
FMN binds, suggesting that the more distant a residue is from FMN, 
the less stabilized it is. However, no significant correlation between 
AAGop the distance a residue has to FMN exists (Supplemen- 
tary Fig. S2). Three peripheral residues become less protected 
against exposure to water on incorporation of FMN (that is, Tyr47, 
Glyl32 and Serl38, coloured dark blue in Fig. 7). Significant FMN- 
induced stabilization of the protein occurs for 54 residues (coloured 
green, yellow, orange and red in Fig. 7), because the corresponding 
AAGop-values exceed 0.5kcalmol~^ (Supplementary Table SI). 

4 



H/D exchange reports cofactor-induced protein stabilization. 

Because of bound FMN, the stability against global unfolding of 
flavodoxin (that is, AG^.uiflavodoxin)) is higher than that of apo- 
flavodoxin. This stability difference depends on the concentration 
of free FMN^^. 

AG^_jy (flavodoxin) = AGjy^_|y (apoflavodoxin) 
.i.Tln[l.»] 

in which R is the gas constant, T the absolute temperature, and 
[FMN] the concentration of free FMN in solution. H/D exchange 
experiments are done using non-reconstituted flavodoxin with a 
concentration of about 2.75 mM and consequently, using a of 
(3.82±0.12)-10"12m, the concentration of FMN not bound to 
protein is 102 nM (that is, binding site saturation is > 99.996%). 
According to equation (1), cofactor binding to apoflavodoxin 
contributes 6.0 kcalmol" ^ to flavodoxins global stability. 

One would expect that the residues with the largest AAGop-values 
reportthe magnitude of cofactor- induced stabilization ofholoprotein. 
Indeed, ten residues of flavodoxin show a stability increase compared 
with apoflavodoxin that equals within error 6.0kcalmol~^ Stabili- 
zation with 5.0 < AAGop < 7.0 kcalmor^ (that is, 6.0kcalmor^ ± 
twice the s.d. of AAGop) happens for Gly8, Ala92, Gly95, Tyrl02, 
Aspl43, Phel46, Vall47 Glyl48, Alal50 and Leul51, which are 
coloured orange and red in Fig. 7. Hence, on binding of FMN, these 
parts of the protein stiffen to such an extent that H/D exchange of 
the corresponding amides from holoprotein hardly happens. How- 
ever, on occasional release of FMN, this exchange does occur from 
apoflavodoxin and is subsequently detected, thereby quantitatively 
revealing the cofactor-induced stabilization of holoprotein against 
global unfolding. Most of these ten amino-acid residues, except for 
Gly8 and Tyrl02, are far from the flavin and experience long-range 
stabilization effects. 

The AAGop -value of one residue of flavodoxin is remarkably high 
(that is, Alal8; AAGop>7.4kcalmori; Supplementary Table SI), 
because it significantly exceeds 6.0 kcal mol~ ^ Apparently, exchange 
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Figure 4 | Typical H/D exchange curves for baclcbone amides of 
flavodoxin. Shown are time-dependent decreases of maximal intensities 
of cross-peaks arising from backbone amides of (a) Arg24 (k^^ = (4.12 
±0.05)-10-4s-''), (b) Ala45 (/Cex = (3.27±0.02)-10-5s-^), (c) Phe49 
aex = (4.14±0.01)-10-6s-^), (d) Val125 aex = (5.16±0.04)-10-^s-^) 
and (e) LeullO (/Cgx^ (2.15 + 0. 03)-10"^s"^), respectively. Notice the 
hugely different time-scales in panels a-e. Data are extracted from ^H-^^N 
HSQC spectra. 



of the backbone amide of Alal8 happens less frequently from apo- 
protein that is temporarily formed on release of FMN from flavo- 
doxin, than occurs from equilibrated apoflavodoxin. We propose 
that Ala 18 is involved in slow structural rearrangements that lead to 
relaxation of apoflavodoxin after release of FMN. 

Intermediate stability increases, with AAGop ranging in between 
0.5 to 5.0kcalmol~^ occur for 43 residues (coloured green and 
yellow in Fig. 7). Because of the stabilizing effects that FMN 
imposes, the corresponding opening processes happen less fre- 
quently than in apoflavodoxin, but more often than release of FMN 
from holoprotein takes place. 

For most residues that are directly involved in intramolecular 
interactions with FMN, no H/D exchange rates could be measured 
in apoflavodoxin (Fig. 6a; Supplementary Table SI) and thus no 
AAGop-values were determined. H/D exchange of these residues 
in apoflavodoxin was not detected, because dynamic exchange 
of the flavin-binding site between different conformations on the 
micro- to milH-second timescale causes many of the corresponding 
NMR cross-peaks to broaden beyond detection^^. Eighteen of these 
residues have AGop-values in holoprotein of less than 5 kcalmol"^ 
(Fig. 6c; Supplementary Table SI), which is significantly less than 
the free energy associated with FMN release from flavodoxin. 




log(/Cint) \og{kJ 



Figure 5 | pD-dependence of hydrogen exchange, (a) Simulated 
pD-dependent exchange curve of an exchangeable proton. The intrinsic 
exchange rate constant, /cj^t, is 1 s"^ at pD 6; the rate constant for local 
opening of the protein, /cQp, is ^x^O~^s~\ and the rate constant for local 
closing of the protein, k^\, is 5s~l An increase of pD by one unit causes 
^jnt to increase tenfold. At low pD-values, exchange rate constant /Cgx 
increases accordingly (area shaded in grey). In this so-called EX2 pD 
regime, the logarithm of depends linearly on pD with a slope of 1. 
Above a certain pD-value, ^jnt is larger than k^\, and exchange enters the 
EX1 regime, in which k^^ curves towards a horizontal line defined by /Cgx 
equalling /Cop. (b) Measured dependence at 25 °C of ^gx o"^ '^int of the 
backbone amides of Leu34 (open circles) and Thr56 (closed circles) 
of flavodoxin, respectively. Leu34 shows typical EX2 behaviour over 
the entire pD-range sampled, because k^^ increases tenfold with 
every increase of pD by one unit. Consequently, AGop can be obtained 
in this pD-range, because AGop= -R7ln(^exAint)- When on increasing 
/cjnt, the value of ^ex curves, as seen here for Thr56, AGop can only be 
determined at low pD-values, and equation (S4) of Supplementary 
Methods can be fitted to the corresponding data (solid line) to 
obtain and k^\. 

as dictated by equation (1). Hence, exchange of these residues 
happens more frequently than release of FMN from holoprotein 
and consequently involves small-scale openings of flavodoxin. 

Rate constants for opening diminish on cofactor binding. 

Cofactor- induced stabiHzation of amino-acid residues against 
exposure to solvent can either be achieved by reducing the opening 
(unfolding) rate of the corresponding protein structure or by accel- 
erating its closing (folding) rate. In case of flavodoxin, 1 1 residues 
are identified for which log(A;ex) has a curved dependency on pD in 
the pD-range used (that is, pD 6.2 to 8.8). For these residues, the 
equation shown in Fig. 5a is used to fit these data (Supplementary 
Methods). This procedure enables determination of rate constants 
for opening and closing of the corresponding residues (Supplemen- 
tary Table S2), with Thr56 as example (Fig. 5b). For five residues 
(that is, Phe6, Ile51, Leu52, Gly53 and Phe94), the rate constants 
for opening and closing are known for apoprotein as well (Supple- 
mentary Table S2). On FMN binding, closing rates reduce 3-to-15 
fold, whereas opening rates of these 5 residues reduce about 5,650- 
fold. Consequently, cofactor binding mainly affects opening events 
of this FMN-binding protein. 

Opening rate of flavodoxin's flavin-binding site. To follow open- 
ing events of the flavin-binding site of flavodoxin, the N(3)H proton 
of FMN (Fig. 8a) is an excellent probe, because its exchange with a 
deuteron of D2O can be measured. To quantitatively interpret H/D 
exchange of N(3)H, one requires knowledge about whether EXl or 
EX2 conditions apply. To distinguish between both situations, we 
determined the intrinsic rate constant for exchange of N(3)H of 
FMN free in solution (that is, A:int(N(3)H)), which was unknown. 
The activation energy of exchange for N(3)H is (13.7±1.1) 
kcalmol"^ and the intrinsic exchange rate of N(3)H at pD 6.3, 
25 °C is (6±2)-103s-i (Methods). 
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Figure 6 | Thermodynamic stabilities of apoflavodoxin and flavodoxin. The cartoon models show AGop at the residue level of (a) apoflavodoxin and 
of (b) flavodoxin, as determined by H/D exchange at pD 6.3, 25 °C. Residues are colour-coded according to the scale bar shown, with AGop-values in 
kcal moTl Residues for which AGop-values could not be determined are shown in grey, (c) Same as (b), but showing only those residues of flavodoxin 
for which AGop is smaller than 5 kcal moTl The cartoon models are generated with PyMOL (Schrodinger, LLC, Palo Alto, CA, USA) using the crystal 
structure of A. wne/and/V flavodoxin (pdb ID lYOB^^), and the FMN cofactor is shown in stick representation. 
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Figure 7 | Binding of FMN thermodynamically stabilizes flavodoxin. 

The cartoon models show the differences in AGop 3t the residue 
level (that is, AAGop) between flavodoxin and apoflavodoxin at 
pD 6.3, 25 °C. The models are rotated by 90° relative to one another. 
Residues are colour-coded according to the scale bar shown, with 
AAGop-values in kcal moPl Residues that are not affected on incorporation 
of FMN (that is, -0.5< AAGop< 0.5 kcal mol"^) are coloured pale blue. 
Residues for which AAGop-values could not be determined are shown 
in grey. 



The exchange rate of the N(3)H proton of FMN in flavodoxin is 
(3.803±0.007)-10-6s-i at pD 6.3 and 25 °C (Fig. 8b). As the cor- 
responding A;int(N(3)H) is large (that is, 6x10^ s~^), we expect that 
N(3)H in flavodoxin exchanges according to the EXl mechanism. 
For EXl conditions to apply, the closing rate A;^ of the microen- 
vironment that surrounds N(3)H in flavodoxin needs to be much 
smaller than A:int(N(3)H). This closing reaction is due to binding 
of FMN to apoflavodoxin and/or due to folding of the microen- 
vironment surrounding N(3)H. Both processes proceed much 
slower than 6x10^ s~^ as the following observations support. First, 
the rate of binding of FMN to apoflavodoxin is smaller than 1 s~ ^ 
(calculated using concentrations of both apoflavodoxin as well as 
free FMN of about 0.1 |aM in a 2.75-mM flavodoxin-containing 
hydrogen exchange sample, and using a rate constant for binding 
of FMN^^ of about 1 jiM" ^). Second, the rate constants for clos- 
ing of flavodoxin (Supplementary Table S2) and of apoflavodoxin^^ 
are well below 6x10^ s~^ Thus, exchange of N(3)H of flavodoxin 
indeed occurs according to the EXl mechanism, and consequently 
at pD 6.3 and 25 °C the opening rate of flavodoxin's flavin-binding 
site is (3.803±0.007)-10-6s-i. 

6 



Discussion 

This paper provides new insights into cofactor-binding interactions. 
We show by using a well-established titration procedure that the 
dissociation constant of flavodoxin is 3.51x10" ^^M (Fig. la). After 
each titration step, part of the free FMN binds rapidly to apofla- 
vodoxin, causing the associated FMN fluorescence to be quenched 
nearly instantaneously to low level. This study reveals for the first 
time that the used titration procedure underestimates the K^-value 
of cofactor-bound protein, because after this initial quenching 
extremely slow conformational events happen. Reconstituted holo- 
protein relaxes with a time constant of ~5 days (Fig. 2b) to a state 
that binds the flavin cofactor about two orders of magnitude better 
than freshly formed flavodoxin. Many residues are involved in this 
relaxation process (Fig. 3). Several backbone amide cross-peaks shift 
with a time constant similar to the one observed by fluorescence, 
whereas some move more slowly. Apparently, the latter relaxation 
phenomena do not affect FMN fluorescence. Relaxed flavodoxin is 
characterized by virtually no fluorescence of bound FMN (Fig. 2a) 
and has a of 3.82x10" ^^M (Fig. lb). Thus, tightening of FMN 
binding leads to an energetically highly favourable protein state 
with picomolar-binding affinity. 

Intrinsic exchange of N(3)H in free FMN is about three orders 
of magnitude faster than intrinsic exchange of backbone amides 
(that is, /:int(N(3)H) equals 6x10^ s'K at pD 6.3 and 25 °C). The 
corresponding activation energy is 1 3.7 kcal mol ~ ^ much less 
than the 17kcalmol~^ required for backbone amide exchange"^^. 
This difference is due to the molecular structure of FMN in 
which N(3)H is chemically attached to two carbonyl groups 
(Fig. 8a), whereas a backbone amide is coupled to only one 
carbonyl. Consequently, above pD 6, exchange of N(3)H of 
protein-bound flavin occurs according to the EXl mechanism 
and is an excellent probe to directly measure opening rates of 
flavin-binding sites. 

In flavodoxin, N(3)H of FMN is tremendously shielded from 
solvent. Compared with free FMN, exchange of N(3)H is slowed 
down by a factor of 1.6x10^ and thus is a rare event. Each opening 
event that makes N(3)H accessible to D2O causes exchange to occur, 
because the EXl mechanism applies. The rate constant of FMN dis- 
sociation from flavodoxin is 1.05xl0~^s~^ (Fig. lb). Hence, H/D 
exchange of N(3)H occurs (3.6 ±0.1) times faster than full release of 
FMN dictates. Consequently, small-scale opening of the microen- 
vironment surrounding N(3)H without dissociation of FMN from 
the protein must occur, and is in fact the main contributor to this 
exchange. 

After nascent apoflavodoxin is released from a ribosome, the 
protein can autonomously fold to its native state^^"^^. During 
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Figure 8 | Hydrogen exchange of N(3)H of f lavodoxin's cofactor. 

(a) Structure of FMN, with N(3)H indicated, (b) H/D exchange curve 
obtained for N(3)H of flavodoxin at pD 6.3, 25 °C. Because of folding 
in the ^^N-dimension of the ^H-^^N HSQC experiment, this intensity is 
negative. The exchange rate is (3.803 + 0. 007)-10"^s"l 



apoflavodoxin folding, no transmission of FMN-binding effects 
happens, because neither unfolded protein nor apoflavodoxins fold- 
ing intermediates interact with FMN^^ Native apoflavodoxin binds 
to FMN firmly. Under reducing conditions, which prevail in the 
cytoplasm, tightening of cofactor binding happens rapidly. 

Because protein- cofactor interactions strongly affect the redox 
potentials of FMN^^~^^, and because the function- altering amino- 
acid residues have been evolutionarily selected to optimize protein 
function, flavodoxin can transfer electrons between various physi- 
ological redox partners. We report here the subtle interplay between 
a protein and its cofactor. By following NMR chemical shifts of 
backbone amides after reconstitution of flavodoxin, we identified 
amino -acid residues that tighten cofactor binding from nanomolar 
to picomolar affinity. By using H/D exchange and NMR spectro- 
scopy, we probed the local free energy differences between apoflavo- 
doxin and fully relaxed flavodoxin and pinpointed amino acids that 
cause picomolar-binding affinity of the cofactor. Both experiments 
show that the flavodoxin sequence has evolved such that cofac- 
tor-induced effects on local structure and stability not only occur 
adjacent to FMN, but also happen in many other regions, includ- 
ing distal residues positioned in the core and at the periphery of 
the protein (Figs 3 and 7). By means of interaction networks, FMN- 
induced stabilization transmits throughout the entire protein. These 
networks within apoflavodoxin and flavodoxin comprise hydrogen 
bonds of a-helices and of the parallel p-sheet, and also involve other 
intramolecular interactions. From a statistical thermodynamic point 
of view, the binding of FMN leads to a redistribution of the native - 
state ensemble. Those states that bind to the cofactor are stabilized 
with respect to those states that do not, causing a change in the 
probability distribution of states^^'^^ This redistribution affects not 
only residues in direct contact with FMN, but also regions linked by 
cooperative interaction networks. Conversely, but fully equivalently, 
these networks of interacting amino-acid residues cooperate to 
firmly bind to FMN. Hydrogen exchange of ten residues, coloured 
orange and red in Fig. 7, reports the FMN-induced global stability 
increase. Thus, residues that stabilize the interaction networks and 
cause tight binding of FMN to flavodoxin have now been experi- 
mentally identified for the first time. Because many more amino 
acids are evolutionarily conserved among flavodoxins, use of align- 
ment of flavodoxin sequences cannot reveal these residues. Two of 
the residues we identified (that is, Gly8 and Tyrl02) are in close con- 
tact with FMN. The other eight amino-acid residues (that is, Ala92, 
Gly95, Aspl43, Phel46, Vall47, Glyl48, AlalSO and LeulSl) are 
far from the flavin. They are dispersed throughout flavodoxin, seem 
irrelevant to its function, and exert long-range stabilization effects 
on flavin binding. 



This study shows that use of NMR spectroscopy and H/D 
exchange can identify amino-acid residues that cause, by means of 
interaction networks, the ultimate tight binding of a ligand. Experi- 
mental identification of these residues not only is required to under- 
stand the thermodynamics of ligand binding, but also is important 
for engineering, designing and evolving ligand-binding proteins. 
We show that the interaction between protein and its ligand is com- 
plex, and that slow relaxation of flavodoxin happens subsequent 
to initial binding of FMN. Similar relaxation phenomena are likely 
relevant for other ligand-binding proteins as well. Ultimately, many 
residues dispersed throughout the flavodoxin structure contribute 
to picomolar-binding affinity of FMN. 

Methods 

Protein purification. The single cysteine at position 69 in wild-type A. vinelandii 
(strain ATCC 478) flavodoxin II was replaced by an alanine to avoid covalent 
dimerization of apoflavodoxin. This protein variant is largely similar to wild- 
type flavodoxin regarding both redox potential of holoprotein and stability of 
apoprotein"^^'"^^. Uniformly ^^N-labelled recombinant C69A A. vinelandii flavo- 
doxin was obtained from transformed Escherichia coli cells and purified"^ This 
protein contains uniformly ^^N-labelled FMN. Subsequently, apoflavodoxin was 
prepared using trichloroacetic acid (TCA) as described'^^''^'^. 

Fluorescence spectroscopy. FMN fluorescence was recorded at 25 °C by using a 
Gary eclipse fluorimeter. FMN was obtained during the TCA-induced preparation 
of apoflavodoxin and further purified by reverse-phase high-performance liquid 
chromatography. During the titration procedure, excitation was at 445 nm with a 
slit of 5 nm and emission was recorded at 525 nm with a slit of 10 nm. During the 
dilution procedure, excitation was set to 450 nm. A sample of 50 nM FMN was 
used to correct for photobleaching and instrumental drift. 

NMR spectroscopy and hydrogen exchange data analysis. Lyophilized protein 
was dissolved in 100 mM potassium pyrophosphate in D2O with 100 |lM 
2,2 dimethyl-2-silapentane-5-sulfonic acid as internal chemical shift reference. 
The final concentrations of protein ranged in between 1.3 to 3 mM. Samples were 
immediately transferred into a 500 MHz NMR instrument. Subsequently, a series 
of HSQC spectra was recorded to detect H/D exchange. The temperature was set 
to 25 °C. The dead time (that is, the difference in time between mixing D2O with 
lyophilized protein and recording of the first HSQC spectrum) was ~5 min. During 
H/D exchange, one-dimensional proton NMR spectra of the sample were also 
acquired. The intensities of the aliphatic resonances in these proton NMR spectra 
do not alter throughout the course of measuring the series of HSQC spectra. This 
observation shows unaltered sensitivity of the NMR spectrometer and integrity 
of the protein sample during our experiments. pD was measured after exchange 
experiments. Spectra were processed and maximal cross-peak intensities were 
analysed as a function of time, according to well-established procedures to 
obtain k^^ (Supplementary Methods). 

Rate and equilibrium constants for FMN dissociation. The increase in flavin 
fluorescence due to FMN release when flavodoxin was diluted in buffer was 
followed at 25 °C. On dilution, the following equilibrium is affected: 

^off 

flavodoxin < ^ apoflavodoxin -I- FMN (2) 
^on 

Return to equilibrium proceeds exponentially as function of time, as it is a 
first-order process, and the corresponding relaxation time informs about /c^ff 
as well as Kd. 

intrinsic rate constant for exchange of N(3)H of FMN. In buffer solution, the 
chemical shift of N(3)H of free FMN can be observed in a ^H NMR spectrum 
when the temperature of the sample is set low, to diminish exchange of this rapidly 
exchanging proton. In addition, one needs a jump-return sequence'^'^ to avoid 
excitation of H2O and subsequent transfer of saturation to N(3)H. At pH 5.22 
and 1.02 °C, N(3)H of free FMN resonates at 11.6p.p.m. 

At pD 6.26 and 25 °C, the conditions at which flavodoxin is probed by 
hydrogen exchange, intrinsic exchange of N(3)H is too rapid to measure directly 
by NMR spectroscopy. However, the corresponding rate constant A;int(N(3)H) 
can be derived indirectly. A selective saturation recovery NMR experimenf*^'"^^ 
was used to determine /cint(N(3)H) at 1.02 °C (verified using 4% methanol in 
d4-methanol) and pH-values of 4.27, 4.72 and 5.22, respectively (Supplementary 
Fig. S3). FMN concentrations were 100 |lM at pH 4.27 and 4.72, and 939 |lM at pH 
5.22, respectively Buffer was 100 mM KPPi, 90% H20/10% D2O. Fitting of a linear 
equation to log(/Cex) versus pH shows that the corresponding slope is 0.93 ±0.1 3, 
which equals unity within error (Supplementary Fig. S3). Consequently, in the 
pH region investigated, exchange of N(3)H is due to catalysis by OH~ . To obtain 
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^int(N(3)H) in D2O, one needs to take into account the difference in ionization 
constants of H2O and D2O, as well as the difference in activities of OH ~ and OD ~ 
(ref 47). In 100 mM KPPi, the activity of OH~ at pH 5.65 equals the activity of 
OD~ at pD 6.26. Thus, at pH 5.65 and 1.02 °C, the intrinsic exchange rate of N(3)H 
is (800 ± 108) s~ ^ To extrapolate A:int(N(3)H) to 25 °C, the activation energy of 
exchange (E^) for N(3)H needs to be known. To obtain this latter value, a selective 
saturation NMR experiment was done at 24.53 °C using a solution of 1 mM FMN 
in 100 mM KPPi pH 4.26. Using these conditions, k^^ is (287.9±22.4) s" ^ How- 
ever, at 1.02 °C, pH 4.27, k^^ is (39.6 ±3.2) s~^ Using equation 3, the activation 
energy of exchange for N(3)H was calculated to be (13. 67± 1.10) kcalmol" ^: 

iCi„,(297.68) = fci„,(274.17K).e-^»'««"297.68)-(l/274.17)) 

with R the gas constant (1.987xlO~3kcalK~ ^ mol" i). 

Using the above results, the intrinsic exchange rate of N(3)H at 25 °C and pH 
5.65 was inferred to be (6±2)10^ s~ ^ As discussed, this rate also applies at pD 6.26, 
25°C, the condition at which the AG^p-values of the backbone amides of flavodoxin 
are determined. 

As mentioned, a selective saturation NMR experiment was used to determine 
the intrinsic exchange of N(3)H. In this experiment, the pulse sequence (relaxa- 
tion delay— saturation pulse— variable delay— observation pulse)n was used to 
determine the apparent longitudinal relaxation rate (i/T^ ^pp) of N(3)H, where n 
is the number of scans acquired. The saturation pulse of 0.15 s was at the frequency 
of N(3)H and was of sufficient power to saturate this proton. A relaxation delay of 
0.5 s was used. During the variable delay period, which ranged from 4 \xs to 170 ms, 
recovery of magnetization of N(3)H occurs. This magnetization was subsequently 
detected by employing an observation pulse and used to derive ^pp- A jump- 
return sequence^^ was used as observation pulse to avoid excitation of H2O. In a 
separate experiment, the H2O resonance was saturated during a period of 2.5 s, 
and subsequently a jump-return sequence was used to detect intensity of the 
N(3)H resonance. In a control experiment, pre-irradiation occurred in an empty 
region of the spectrum, far removed from both water and N(3)H resonances, to 
determine intensity Mq of the N(3)H resonance. Subsequently, the intrinsic rate 
constant for exchange (1/t) of N(3)H was calculated using^^: 

l/T = (l/Ti,,pp)-(l-M3/Mo) (4) 

Ms turned out to be zero under all experimental circumstances investigated by us 
and consequently: 

l/'Z^^l/Ti.app (5) 
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